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Abstract 

This letter presents a remarkably simple approach to the first-principles determination of the 
ergodic-non-ergodic transition in monodisperse colloidal suspensions. It consists of an equation 
for the long-time asymptotic value 7 of the mean squared displacement of the colloidal particles, 
whose finite real solutions signal the non-ergodic state, and determines the non-ergodic parameter 
f{k). We illustrate its concrete application to three simple model colloidal systems, namely, hard- 
spheres, hard-spheres plus repulsive (screened Coulomb) Yukawa interaction, and hard-sphere plus 
attractive Yukawa tail. The results indicate that this is quite a competitive theory, similar in spirit 
to, but conceptually independent from, the well-known mode coupling theory. 

Keywords: colloid dynamics, glass transition, dynamic arrest. 

Esta carta presenta un metodo notablemente simple para la determinacion, de primeros principios, 
de la transicion ergodico-no-ergodico en suspensiones coloidales monodispersas. Dicho metodo 
consiste en una ecuacion para el valor asintotico a tiempos largos, 7, del desplazamiento cuadratico 
medio de las particulas coloidales, cuyas soluciones reales finitas son sinonimo de no-ergodicidad, 
y determinan el parametro no ergodico f{k). Ilustramos su aplicacion concreta en tres modelos 
simples de sistemas coloidales, a saber, esferas duras, esferas duras con interaccion de Yukawa 
repulsiva (Coulombica apantallada) , y esferas duras con interaccion de Yukawa atractiva (fuerzas 
de deplecion). Los resultados indican que esta es una teoria muy competitiva, similar en espiritu, 
pero conceptualmente diferente, a la bien conocida teoria de acoplamiento de modos. 

Descriptores: dinamica coloidal, transicion vi'trea, arresto dinmico. 
PACS numbers: 64.70.Pf, 61.20.Gy, 47.57. J- 
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The most basic and elementary information on the thermodynamic properties of a ma- 
terial is its phase diagram. The description of the gas-liquid transition provided by the van 
der Waals equation of state is the earliest and most paradigmatic example of the construc- 
tion of a phase diagram starting from molecular considerations [1]. The most outstanding 
achievement of statistical mechanics, however, has been the establishing of the microscopic 
version of the second law of thermodynamics, which provides the basis for the systematic 
calculation, given the intermolecular forces, of phase diagrams, by the simple conceptual 
procedure of identifying the phase with the lowest free energy. This very fundamental prin- 
ciple allowed the development of the wide range of methods, approaches, techniques and 
applications of equilibrium statistical thermodynamics 2[]. 

At the same time, one of the most notorious limitations of statistical mechanics has been 
its inability to identify an equally general and simple principle that allows us to describe 
non-equilibrium states of matter, given the molecular interactions. From a practical and 
fundamental perspective, this is quite disturbing, given the fact that a large variety of the 
materials with which we actually interact in ordinary life are not in their thermodynamic 
equilibrium state. Thus, it is permanently important to search for first-principles approaches 
to describe the most elementary properties of such "phases" and the transitions between 
them, even in the context of specific classes of non-equilibrium states. Perhaps the simplest 
of them refer to the states that result when kinetic barriers prevent a material from reaching 
its thermodynamically most stable state, thus being trapped in dynamically arrested states, 
such as glasses or gels. One then would like to have first-principles criteria to predict the 
location of the boundary between the ergodic fluid phase and such arrested non-ergodic 
states. Contrary to its thermodynamic analog, in this case we only have a single theory of 
this sort, namely, the mode coupling theory (MCT) of the ideal glass transition 0, 4], some 
of whose predictions have found beautiful experimental confirmation . Unfortunately, 
its criterion to decide in a practical manner if a system is in an ergodic or in an arrested state 
is somewhat obscured by the conceptual complexity of this theory and by the many other 
issues addressed js]. Thus, the need exists of simpler and more straightforward approaches 
that focus on this specific and important issue. The main purpose of this communication is 
to present a new first-principles prescription to locate the ergodic-non-ergodic transition in 
monodisperse colloidal dispersions in a remarkably practical and simple manner, completely 
independent of the MCT and its recent variants p, Isj. 
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FIG. 1: Real solutions 7 of Eq. ([I]). Below 



0.563, this equation has no real solutions. Above 



two solutions appear, illustrated by the two branches that bifurcate at (pg. The branch for 



which 7 decreases with (f) (solid line) corresponds to the physical solution of the glass state. In 
the inset we plot the functional $ [7; S] as a function of 7 for = 0.50, 0.55, 0.60, and 0.65 (from 
bottom to top). 

This prescription consists of an equation for the most elementary "order parameter" 
associated to the transition from an equihbrium fluid phase to a glass or gel state. This is 
the long-time asymptotic value, 7, of the mean squared displacement (msd) of the particles in 
the suspension. In the arrested states, this parameter is finite, representing the localization 
of the particles, whereas in the ergodic phase, it diverges. Such an equation reads 



[S{k)-l]'X\k)^ 



dkk' 



[X{k)S{k) + k^-l] [A(A;) + F7] 



(1) 



where S{k) is the static structure factor of the system, \{k) = [1 + (/c/fcc)^]"^, with k^. being 
the position of the first minimum of S{k), and n is the particle number concentration. The 
very form of this criterion exhibits its simplicity: given the effective inter-particle forces, 
statistical thermodynamic methods allows one to determine S{k)^ and the absence or ex- 
istence of real solutions of this equation will indicate if the system remains in the ergodic 
phase or not. The next objective of this letter is to illustrate the practical application of 
this criterion, in the context of three simple model colloidal systems, and to explain very 
briefly its conceptual origin. 

To illustrate the use of this criterium, let us notice that the left-hand side of Eq. ([T]) 
is a functional of S{k) and an ordinary function of 7, which we denote by $ [7; S]. Thus, 
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for a fixed state [i.e., fixed S{k)], this equation may be solved by plotting $ [7; S*] as a 
function of 7, to see if it crosses unity, and for which value(s) of 7 it does so; notice that 
this functional must vanish in the limits of large and small 7. This very simple procedure 
leads to the determination of the finite solutions of Eq. ([T]), as illustrated in Fig. [T]for the 
hard sphere (HS) system with the static structure factor S{k) given by the Percus-Yevick 
(PY) approximation with the Verlet-Weis (VW) correction [10]. The inset of Fig. [1] 
exhibits the dependence of the functional on 7 for various volume fractions. Clearly, below 
a threshold value (pg, $ [7; S] remains below 1.0 for all 7, and hence, there are no real 
solutions; thus, the system must be in the ergodic state for < 0g. Above (j)g there are two 
real solutions, the smallest one corresponding to the glass state, since in the glass the msd 
must decrease with n. In this manner, we determine this threshold value to be (pg = 0.563. 
This quantitative prediction is closer to the experimental values than those of the original 
MCT or its extensions. 

of the self- 
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13j. This 



The criterion above emerges from the long-time asymptotic analysis [11 
consistent generalized Langevin equation (SCGLE) theory of colloid dynamics 
theory was originally developed to describe the dynamic properties of colloidal dispersions. 
Thus, it allows the calculation of properties such as the msd of the particles, or the interme- 



diate scattering function F{k,t) and its se// component, Fs{k,t) (ij, ll5|, given the effective 



pair potential u{r) between colloidal particles. From the same long-time analysis one can 
also derive [ll| a remarkably simple expression for another experimentally important prop- 
erty of the glass states, namely, the non-ergodic parameter f{k) = F{k,t 00)/ S{k). Such 
an expression reads 

_ Xik)Sik) 

X{k)S{k) + k^Y ^ ^ 

with 7 being the physical solution of Eq. ([1]). Let us notice that this result only involves 
S{k) and X{k) as static inputs. We may illustrate its use by applying it to the HS system 
above, for which we just found (pg = 0.563. Right at (pg, there is only a single solution for 
7, namely, 7 = 1.06 x 10~^cr^, with a being the HS diameter. This solution for 7 may then 
be employed in Eq. ([2]) to determine the non-ergodic parameter f{k). Fig. 2 compares this 
prediction for f{k) with the experimental data of Ref. 5|. A similar comparison is also 
included in this figure, corresponding to a higher volume fraction, (p = 0.58, and to the 
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FIG. 2: Non-ergodic parameter f{k) (solid lines) calculated with Eqs. ([T]) and ([2]) for the HS system, 
with S{k) given by the PY approximation with the VW correction, at the theoretically-predicted 
ideal glass transition, (pg = 0.563 (lower curve), and at the larger volume fraction (p = 0.58 (upper 
curve; shifted upwards one unit, for visual clarity). Symbols are the experimental data of Ref. 
(filled circles) and Ref. [6|] (empty circles). 

experimental data of Ref. 6|]. 

The next example refers to dispersions of charged colloidal particles, with effective pair 
potential modeled by the HS plus the repulsive screened Coulomb potential, Pu{r) = 
K exp[—z{r /a — 1)]/ (r/cr) for r > a. The interaction parameters z and K are, respectively, 
the inverse Debye screening length (in units of a), and the intensity of the pair potential 
at hard-sphere contact (in units of the thermal energy (3~^ = ksT). We employed as the 
static input of Eqs. ([1]) and ([2]) the experimentally-determined static structure factor of a 
sample with a = 272nm and (j) = 0.27, provided in Ref. Q]. For this, we used the hyper- 
netted chain approximation [2] to provide a smooth fit of the data, leading to 2 = 3.14 and 
K = 11.56. The solution of Eq. ([I]) is 7 = 0.00293o-2, and the results for f{k) from Eq. ([2]) 
are compared with the corresponding experimental data in Fig. 3. Given that no fitting 
parameters are involved in the theoretically-predicted non-ergodic parameter f{k) in Figs. 2 
and 3, we consider that the overall quantitative accuracy of our theory is quite remarkable, 
and is certainly better than that of the MCT. Of course, one can easily determine the full 
liquid-glass "phase" diagram in the space {K, z, (j)), but this is not discussed here. 

The last illustrative example refers to another simple model of a colloidal dispersion, 
this time involving the hard-sphere plus an attractive Yukawa (HSAY) potential, /5u(r) = 
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FIG. 3: Theoretical prediction of the non-ergodic parameter f{k) (solid line) calculated with Eqs. 
([1]) and ^ for the hard-sphere plus the screened Coulomb tail, with (j) = 0.27, z = 3.14, and 
K = 11.56. The symbols are the experimental data of Ref. [gI 
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FIG. 4: Theoretical glass transition line in the plane (T*, 0) for the HSAY system, calculated with 
Eq. ([I]) (thick solid and dashed curve), and from the MCT (symbols; taken from Ref. 

Q)- The 

lighter (solid and dashed) line corresponds to the freezing line, according to the HV criterium. 
The dotted curve corresponds to the spinodal line. The inset displays the same information in the 
{K, (p) plane 



—Kexp[—z(r/a — l)]/(r/o"), modeling depletion forces [16|]. In this case we only emphasize 
some qualitative aspects of the predictions of our theory, obtained using the mean spherical 
approximation (MSA) [itI for S{k). Fig. 4 presents the ergodic-non-ergodic transition line 
of this system for 2; = 20 in the plane {T*,(f)) of state parameters, with T* = being 
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a reduced temperature. For reference we also plot the line where the maximum of S{k) 
reaches 2.85 (the freezing line, according to the Hansen- Verlet (HV) criterium [3]), as well 
as the spinodal curve. To illustrate the qualitative similarity between our results and those 



of the MCT, we have also included the MCT transition line [19|. The inset re-plots the 
same curves, now in the plane {K,(f)), to illustrate the experimentally observed [l6[] shape 
of the glass transition line at its high density end. This describes the re-entrant behavior 
of a system upon lowering its temperature at fixed volume fraction, from a hard-sphere 
glass to an ergodic state, followed by the re-entrance from the ergodic state to a second 
("attractive") glass state. 

This completes the illustration of the practical use of the general results in Eqs. ([T]) and 
([2]). Let us now comment on their physical origin. Of course, the physics behind them 
is the physics behind the full self-consistent theory from which they derive. Hence, let us 
summarize the four distinct fundamental elements of this theory. The first consists of general 
memory-function expressions for F{k,t) and Fs{k,t) derived with the generalized Langevin 



equation (GLE) formalism UJ, ll2|, which in Laplace space read 



^ + l+C{k,z) 



Psik,z) = -^^^, (4) 



1+Cs{k,z) 

where Dq is the free- diffusion coefficient, and C{k,z) and Cs{k,z) are the corresponding 
memory functions. 

The second element is an approximate relation between collective and self-dynamics. We 
approximate the difference [C{k,t) — Cs{k,t)] by its exact shor t-t ime / large- k limit, thus 



defining what we refer to as the "additive" Vineyard- like approximation 



12|, 



C{k,t) = Cs{k,t) + [C^^-^^(A;,t) - C|^^^(fc,t)]. (5) 
In this equation, C^^^^{k, t) and C§^^^{k, t) are the exact short-time expressions for these 



memory functions, which also define the so-called 15|, l20[ single exponential (SEXP) ap 



proximation, and for which well-established expressions, in terms of equilibrium structural 



properties, are available 
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The third ingredient consists of the independent approximate determination of Fs{k,t) 
[or Cs{k,t)]. One intuitively expects that these /c-dependent self-diffusion properties should 
be simply related to the properties that describe the Brownian motion of individual parti- 
cles, just like in the Gaussian approximation [14], which expresses Fs{k,t) in terms of the 
mean-squared displacement W(t) as Fs{k,t) = exp[—k'^W(t)]. We introduce an analogous 
approximate connection, but at the level of their respective memory functions. The mem- 
ory function of W{t) is the so-called time-dependent friction function A({t). This function, 
normalized by the solvent friction (q, is the large wave-vector limit of Cs{k,t). Thus, we 
interpolate Cs{k,t) between its two exact limits, namely. 



Cs{k,t) = Cl^''''{k,t) + [AC{t) - A(A;), (6) 



where A(*{t) = A({t)/(Q. The fourth ingredient o: 
this property, also derived with the GLE approach 



our theory is a general expression for 



2l|, namely. 



k[S{k) - 1] 



1 2 



Fik,t)Fsik,t). (7) 



S{k) 

Eqs. (l3])-(I7j) constitute the SCGLE theory of colloid dynamics. Besides the unknown 
dynamic properties, it only involves the static structural properties S{k), C'^^^^{k,t) and 
Cg^^^{k,t), determined by the methods of equilibrium statistical thermodynamics. Con- 
cerning the interpolating function X{k), phenomenological arguments were given that 
led to the definition given above [immediately after Eq. ([1])]. Although no fundamental basis 
is available for this choice of X{k), this definition is universal (in the sense that it is the same 
for any system or state), and renders the resulting self-consistent scheme free from any form 
of adjustable parameters. 

The results in Eqs. ([1]) and ([2]) were derived from the long-time analysis of the SCGLE 
scheme in Eqs. ©-([T]). We have, however, also solved numerically the full self-consistent 
theory for the systems discussed here. Of course, we confirmed, with this lengthier method, 
the quantitative results obtained in a much more economical manner from Eqs. ([T]) and ([2]). 
The numerical solution, however, provides the whole scenario of the relaxation processes. 
On the basis of such results, which we shall discuss separately, we may say that the general 
scenario of the glass transition provided by the present theory is consistent with the available 
experimental data, and is quite similar to that provided by the MCT. We should point out. 
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however, that for the HS system, our theory does not have to appeal to any sort of re-scahng 
of the volume fraction, as the MCT is forced to invoke due to its numerically low predicted 
glass transition volume fraction (^^ = 0.52). In fact, none of our theoretical results presented 
here involved any form of adjustable parameter. Thus, we conclude that in many respects, 
the present theory of dynamic arrest is a sound and competitive theoretical description of 
dynamic arrest in colloidal systems. 
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